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Abstract 

We propose analytic approximations of chiral SU (3) amplitudes for the extrapolation of 
lattice data to the physical meson masses. The method allows the determination of NNLO 
low-energy constants in a controllable fashion. We test the approach with recent lattice 
data for the ratio Fk/F-^ of meson decay constants. 
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1. Introduction 

In recent years, lattice QCD has made enormous progress in the hght quark sector (see, e.g., 
Ref. [1]). State-of-the-art lattice studies employ quark masses corresponding to pion masses 
as low as 200 MeV, with kaon masses close to the physical value. 

Extrapolation to the physical meson masses is performed in different ways. On one side 
of the theory spectrum, there are various smooth extrapolation formulas with more or less 
theoretical motivation. On the other side, the most sophisticated extrapolations are based on 
chiral perturbation theory (CHPT), the effective field theory of the standard model at low 
energies [2l|3]. As the meson masses continue to approach the physical values in future high- 
statistics simulation^, even simple-minded polynomial approximations will allow predicting 
physical quantities with ever better precision. However, a lot of information about QCD is lost 
in this way. On the other hand, CHPT provides the correct analytic structure of amplitudes in 
terms of several a priori undetermined constants, the so-called low-energy constants (LECs), 
which are independent of the light quark masses by definition. 

Many of the higher-order LECs are difficult if not impossible to extract from actual ex- 
perimental data. Lattice calculations offer a new environment for determining LECs because, 
unlike nature, the lattice physicist can tune the quark masses. For the chiral practitioner, it 
is then an advantage rather than a drawback that present lattice studies work with different 
meson masses larger than the physical values. 

Many lattice groups use next-to-leading-order (NLO) CHPT results for the chiral extra- 
polations. As by-products, several LECs at this order, 0{p^), have actually been determined 
this way (see, e.g., Ref. [5]). On the other hand, state-of-the-art NNLO CHPT results have 
only very recently been used for the interpretation of lattice data [6HH]- There are good 
reasons why lattice physicists have generally ignored available NNLO calculations so far: the 
results are quite involved and, what is even worse, they are mostly available in numerical 
form only, at least for chiral SU(3) (for a review of NNLO results, see Ref. [9]). 

We propose in this note analytic approximations for NNLO CHPT amplitudes for chiral 
SU(3) that are more sophisticated than the double-log approximation [lOj, yet much simpler 
than the full numerical expressions. We first recapitulate CHPT to 0{p^) for the generating 
functional of Green functions [11]. The form of this functional suggests analytic approxima- 
tions of p^-amplitudes that are scale independent in contrast to the double-log approximation. 
Moreover, they include all leading and next-to-leading contributions at large Nc- In a first 
exploratory study, we then apply our approximation to the very recent results of the BMW 
collaboration [12] for the ratio Fk/F^^. 

2. Chiral perturbation theory to 0{p^) 

The most compact representation of CHPT in the meson sector is in terms of the generating 
functional of Green functions Z[j] [2113] . In the sequel, we suppress the dependence on external 
fields j. Analogous to the chiral Lagrangian, the generating functional permits a systematic 
chiral expansion: 

Z = Z2 + Z4 + Zg + ... (1) 

The NNLO functional Zq of 0{p^) is itself a sum of different contributions shown pictori- 
ally in Fig. [TJ In the following, we recapitulate and reformulate the general treatment of 

^ A very recent simulation [J already uses physical quark masses. 
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Figure 1: Skeleton diagrams for the generating functional of 0{p^). Simple dots, crossed 
circles, black box denote vertices from leading-order, NLO, NNLO Lagrangians, respectively. 
Propagators and vertices carry the full tree structure associated with the lowest-order La- 
grangian. 



renormalization at 0{p^) [TT]. In addition to tree diagrams of 0{p^) (diagram g), there are 
two classes of contributions requiring separate treatments: irreducible (diagrams a,b,d) and 
reducible (diagrams c,e,f) contributions. 

With dimensional regularization, the irreducible diagrams have both double- and single- 
pole divergences. Moreover, the single-pole divergences of each irreducible diagram are in 
general non-local. Renormalization theory guarantees, however, that the sum of the three 
diagrams has only local divergences |1HI13]. i.e., polynomials in momenta and masses in 
momentum space. Chiral symmetry guarantees that these divergences can be absorbed by 
the LECs of 0{p^) via diagram g. In this process an arbitrary renormalization scale is 
generated. The sum of diagrams a,b,d,g is then finite and can be written in the form (details 
of the derivation will be given elsewhere 



a+b+d+g 



4r«L-r(2)L2 + 2r(^)(/.)L 



Oa 



+ 



(4vr) 



Hi{x-M) + 



(4vr)^ 



.K{x-M) 



(2) 



The structure ^ holds for chiral SU (n) in general but we have already used the notation 
for re = 3. The monomials Oa{x) {a = 1, . . . ,94) define the chiral Lagrangian of 0{p^) [15j 
with associated renormalized LECs C^di), the L1{(£) {i = 1, . . . , 10) are renormalized LECs 
of 0{p^) with associated beta functions Fj [3] and the coefficients Fi^\ fI^-* and fI^^ are listed 
in Ref. [TT]. Repeated indices are to be summed over. Fq is the meson decay constant in the 
chiral SU (3) limit. The chiral log 

L = ^lnMV/.^ (3) 

involves an additional (arbitrary) scale M but ^s well as the total generating 

functional are independent of both fi and M. Hi{x;M) are one-loop functionals associated 
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with diagram d whereas the two-loop contributions (except for the chiral logs) are contained 
in the functional K{x;M). Unlike Oa{x), the functionals Hi{x;M), K{x]M) are non-local. 
The scale independence of the functional ([2]) can be derived with the help of renormalization 
group equations for the LECs CKn) [11] and L\{ijl) [3], using also relations between the 



coefficients ri^^ and ri^^ 

As shown in Ref. [11], the sum of reducible diagrams c, e, f gives rise to a finite and scale 
independent functional with the conventional choice of chiral Lagrangians. It can be written 
in the form 



7C+e+f 



- ^L] Pi^^{x) + F^{x;M) 



d'^x d^y 



Ga,f3{x,y) 



(4) 



Indices a, (3 run from 1, . . . , 8 (octet of pseudoscalar mesons), the Pi^a{x) are local functionals, 
the one-loop contributions of diagrams c,e are contained in the non-local functionals Fa{x] M) 
and the Ga,/3{x,y) are (functional) propagators. 

The complete generating functional of O(p^) is then given by the sum 



^a+b+d+g , /7c+c+f 
^6 + -^6 



(5) 



Once again, it is independent of both scales fi and M. 



3. Analytic approximation for chiral SU{3) 



As emphasized in the introduction, the genuine two-loop contributions contained in the func- 
tional K{x;M) are usually only available in numerical form for chiral SU{3). On the other 
hand, the one-loop contributions can be given in analytic form and the dependence on meson 
masses is manifest. For the chiral extrapolation of lattice results, we therefore suggest the 
following approximate form for the functional of 0{p^): 



+ 



d\ 



1 



+ 



(4vr) 



Clin) + 4r« L - ri^) l2 + 2r(^)(^)L 



Oa{x) 



Hi{x;M) 



d'^x d^y 



+ 2 L[(^) 



-^L]PiAx) G^Ax,y)Fp{y-M) 



(6) 



In contrast to the generalized double-log approximation [10], which can be recovered by 
setting the coefficients ri^^ and the functionals Hi{x;M), Fa{x; M) to zero, the analytic 
approximation ^ is scale independent. This is an important asset for a reliable determination 
of renormalized LECs. 

In addition to tabulated quantities [11], the finite and scale independent one- loop func- 
tionals Hi(x; M) (i = 1, . . . , 10), Fa{x; M) (a = 1, . . . , 8) must be determined from diagrams 
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d and e. If available, the corresponding amplitudes can be read off from existing calculations 
by collecting all amplitudes linear in the LECs L\{^). 

Another attractive feature of ^ is its large- A'^c behaviour. It comprises all leading (gener- 
ically: Ca, LiLj) and next-to-leading contributions (Lj x 1-loop). Of the NNLO terms, it 
contains at least all chiral logs. 

The amplitude for a given observable corresponding to the scale invariant functional ([6]) 
can now be determined in four steps. 

1. Calculate all tree- and one-loop diagrams, i.e., the contributions from diagrams d,e,f,g 
in Fig. [TJ In many cases of interest, these amplitudes are already available in the 
literature [9]. 

2. In the tree-level amplitude of 0{p^) (diagram g), replace the LECs C^(/i) by 

c:{f^) c:{^i) + ^ (4r« L - r(2) + 2rW(/x)L) . (7) 

We recall that the combination on the right-hand side of ([7]) is scale invariant. 

3. Collect all contributions linear and bilinear in the LECs L^{p,) in the remaining am- 
plitude (diagrams d,e,f) and extract the chiral logs. The products L'[{^) L from the 

irreducible parts must match the terms t[^\^)L in ([7]). After performing this check, 
set all chiral logs L = in this subset of terms (diagrams d,e,f), which amounts to 
replacing by M in the one-loop functions. 

4. Replace the bilinears L\[ii)U-[^) (due to reducible contributions: diagram f) by the 
scale invariant expressions 

^{^,)^{^,) (l[(/.) - l^mi^) - ^l) . (s) 

Finally, in the remaining terms linear in the L^fj,) (originating from diagrams d,e) 
replace 

LK/x) - . (9) 

The resulting scale invariant amplitude corresponds to the functional ([6]) . The approximation 
consists in dropping K(x;M) and the terms bilinear in Fa{x;M) in the exact functionals 
dl]) and dH), introducing a dependence on the scale M. This scale parametrizes the two- loop 
contributions not contained in ([6]). Transforming the one-loop functionals Hi{x; M),Fo,{x; M) 
back to Hi{x; fi), Fa{x; fi), the only M-dependence resides in the chiral logs. The remaining 
(single and double) chiral logs can then only be due to the two-loop contributions because all 
other contributions are correctly included in the approximate functional ([6]) and are therefore 
independent of M. Experience with the double-log approximation [TO] suggests that M is 
naturally of the order of the kaon mass in SU (3) calculations. 

While the above approximation is motivated by large Nc, some of the terms not included 
in the approximate functional ([6]) have a relatively simple analytic form (products of one- loop 
amplitudes from diagrams a,c in Fig. [1]). In practice, inclusion of those terms may improve 
the accuracy of the approximation for certain observables. We will come back to this issue in 
Ref. [H]. 
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4. Application to lattice data for Fk/F.,^ 

We apply the analytic approximation ([6]) to the ratio Fk/F^^ of meson decay constants. 
Fk/F-k is well suited for an exploratory study for at least two reasons. 

• At the scale // = 0.77 GeV, the genuine two-loop contribution amounts to — 0.005 for 
physical masses, i.e., half a percent only [6l[T6]. 

• The detailed results of the BMW collaboration [12j provide an ideal laboratory for 
testing our approximation. 

The approximate form of Fk jF^^ is given in the Appendix where all masses are lowest-order 
masses of 0{p^) [16yi7|. Since we work to 0{p^) the masses in must be expressed in terms 
of the lattice masses to O(p^) [3]. The chiral limit value Fq is deduced from the experimental 
value F^ = 92.2 MeV and physical meson masses, using again the relation to 0{p^) [3J. 

Here we are mainly interested in getting information on the LECs of 0{p^). Only two 
combinations of LECs appear: C14 -|- C15 and C15 -|- 2 C17. Most of the Li also contribute to 
Fk/ F-K- Several determinations of the L[ are available in the literature [ISpiQj. All fits of the 
L[ to 0{p^) need to make some assumptions about the C^, in particular for extracting L5 
from Fk/Ft^. L5 is the only LEG contributing to Fk/Ft^ at 0{p^). Moreover, L5 appears at 
Oijf') to leading order in l/Nc- This suggests to fit the lattice data with the three parameters 
Lg, -|- and C^^ + 2C'[j. For the remaining LECs of 0{p^) we adopt the values of fit 
10 of Ref. [18]. Since our approximation is scale independent we may choose the conventional 
scale IX = 0.77 GeV. 

Restricting the data sample of the BMW collaboration |12j to simulation points with 
< 450 MeV, we are left with 13 data points. For this exploratory study, we take only 
the statistical errors of Fk/F^^ into account. After all, our main purpose is to investigate the 
capacity of lattice data for the determination of LECs but not to compete with the detailed 
analysis of Ref. |12] . 

Fitting the 13 data points with Eq. (fT3|) . we obtain for the LECs at = 0.77 GeV 

(0.76 ± 0.09) • 10^3 
(0.37 ± 0.08) • 10"^ GeV"2 

(1.29 ± 0.16) • 10~^ GeV"2 . (10) 
The three parameters are strongly correlated, with correlation coefficients indicated below. 





LI 
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CI, + 2 CI, 


0.69 


-0.87 


LI 




-0.95 



Taking these correlations into account, Fk/F.,^ for physical meson masses is found to be 

Fk/F^ = 1.198 ± 0.005 , (12) 

comparing well with the result Fk/F-,^ = 1.192(7)stat(6)syst of Ref. [12]. We stress once more 
that our errors take only the statistical errors of the lattice values for FkIF.^ into account. 
The same word of caution applies to x^/dof = 1.3 for the quality of fit. 
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Although the functional ([6]) and therefore Fk /F^^ in (jl3p are independent of the renormal- 
ization scale // there is a residual dependence on the scale M of the chiral logs. As announced 
before, we adopt the natural choice M = (lattice value). Varying this scale by it 20 %, 
both Fk/Ft^ and remain practically unchanged while the LECs of 0{p^) vary within two 
standard deviations. Note that this range for M includes M^, which is given by the Gell- 
Mann-Okubo mass formula in the contribution of 0{p^) and is therefore always less than 1.2 
Mk for the meson masses under consideration. Taking M too low would enhance the chiral 
logs too much for an SU (3) observable. The sensitivity to the scale M could be substantially 
reduced if lattice simulations would use strange quark masses lighter than the physical value. 
In such a scenario the convergence properties of chiral SU{3) could be improved altogether. 

For the case at hand, Fk/F.,^ and are insensitive to the approximation made, with 
uncertainties determined by lattice errors. The situation is opposite for the LECs of 0{p^): 
here the approximation errors definitely exceed the lattice errors. The dependence on the 
other LECs of 0{p^) must be taken into account in addition. 

Since C15 is subleading in 1/Nc our fit determines essentially C14 and C17 [6]. Although 
the values depend of course on the input for the L[ we have found generically that both Cl^ 
and are positive and smaller than 10~^ GeV~^, always taken at the usual scale fi = 0.77 
GeV. Comparing with resonance exchange predictions [20], our results indicate that multi- 
scalar exchange is important for these LECs. With single resonance exchange only, we would 
have = Cyj < instead [20]. Our result for Lg lies in the range covered by other NNLO 

fits [Ullg]. 

We cannot compare directly with the results of Bernard and Passemar [6j for Cl^ + Cl^ 
and CJg + 2C[7. First of all, the results of the BMW collaboration were not yet available for 
their analysis and, what is probably more important, the value of was taken as input in 
Ref. [6]. For the reasons given earlier and because of the strong (anti-)correlations found we 
consider it more appropriate to fit L5 together with the LECs of 0{p^). Generically, we find 
somewhat bigger values for C^^ + Cl^ and Cl^ + 2C[7 than in Ref. [6] . 

5. Conclusions 

Starting from the structure of the generating functional of Green functions to O(p^), we have 
proposed analytic approximations for chiral SU{3) amplitudes that require the calculation of 
tree-level and one-loop diagrams only. The result serves two purposes: 

• It provides flexible and user-friendly extrapolation formulas for lattice data. 

• It allows for the determination of higher-order LECs that are otherwise difficult to 
extract from experimental data. 

The approximate amplitudes are independent of the renormalization scale, a prerequisite for 
a reliable determination of LECs. In addition to including all chiral logs, the amplitudes 
contain all leading and next-to-leading terms in the expansion. 

The approach will therefore be especially useful in cases where the genuine two-loop 
contributions are small, compatible with the large- A'^c counting. The ratio Fk/F-,^ is an 
interesting observable with this property. Fitting the approximate expression for Fk/F-,^ to 
recent lattice data, we obtain a value for Fk/F-,^ in agreement with the detailed analysis 
of Ref. p^j. Both Fk/Ft^ and are insensitive to the approximation made. The LECs of 
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0{p^), C[4 + C[5 and Cl^ + 2Clj, are consistent with expectations but subject to uncertainties 
exceeding the lattice errors. 

Although the present study is mainly of exploratory nature we consider the results signif- 
icant enough to warrant further investigations along these lines [14] . 
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Appendix: Approximate result for Fk/Ft^ to 0{p^) 

In the approximation defined by the functional ([6]), Fk/F^^ assumes the following form: 

Fk/F^ = l + R^ + R^ (13) 

F^R^ = A{Ml-M^)L5-5A{M^,fi)/8 + A{MK,f-i)/4 + 3A{M^,fi)/8 

F^Re = 8Fo'(^i-^')(2Mi(Ci4 + Ci5) + M2(Ci5+2Ci7)) 

+ (Mi - m2) (-32 (m2 + 2 M|.)L4 L5 - 8 (3 + mI)lI 

+ (25M2 + 17M|.)LV32) 

(Mi-_M2) 
(47r)2 



+ f._,2 (-2 (M^' + Ml)L, - {Ml + Mi)L2 - (5 + M^Lg/lS 



+ 6 (M^ + 2 M|:)L4 + (14 + 22 mI)L^/3 - 12 {Ml + 2 M|:)L6 
+ lQ{Ml - mI)Lt -A{Ml + bMl)L^ + (313 + 271 M|-)L/288) 
+ 5:4(M,, ^)V8 - A{Mk, ^Ji?/S + A{M^,iJi) A(Mi^, ^)/16 
- 3l(M,, ^) A(M^, /i)/8 - 3 A(Mx, ;u) A(M^, ^)/16 
+ :4(M^,/i) {AMlLi + lOM^Ls + I3M2L3/2 + 10 (M^ + 2M|:)L4 
+ (19 Ml-b Ml)Lrj2 - 10 (M^ + 2 M^Lq - 10 M^Lg 

- (361 + 131 Mi)L/288) 

+ :4(Mi^ , /i) (-4 mILi - 10 M|:L2 - 5 M^L^ - 4 (M^ + 2 M|:)L4 

- {Ml + Ml)L^ + 4 (M2 + 2 Mi)L6 + 4 M^Lg + (59 Ml + 115 M|.)L/144) 
+ A{M^,fi){Ml - Ml)/M^ {-9MlL7 -SMILs + 5MIl/32) 

+ A{M^,f,) {{Ml/2 -2Ml)L^ -Q{mI + 2MI)U - {7 MI + 23MI)L^/Q 
+ 6 (M2 + 2 M|.)L6 + 3 (3 MlM^/M^ - 7 + 4 M^Ly 
+ 3 {MlM^/M^ -3Ml + i Ml)Ls 
-(15 M^Mi/M^ - 44M2 - 19Mi)L/96) 

The abbreviations = Ll{fj.), Ca = C]^{lJ') have been used for a compact representation. 
The masses are the lowest-order meson masses of 0{p'^), Fq is the meson decay constant in 
the chiral SU{3) limit and the chiral log L is defined in ([3|). The loop function A{Ma,fJ-) is 
defined as 

3(M.,M)=(^logi|. (14) 
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